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ABSTRACT 

We have analysed the long- and short-term time dependence of the pulse arrival times 
and the pulse detection rates for eight Rotating Radio Transient (RRAT) sources from 
the Parkes Multi-beam Pulsar Survey (PMPS). We find significant periodicities in the 
individual pulse arrival times from six RRATs. These periodicities range from ~30 
minutes to 2100 days and from one to 16 independent (i.e. non-harmonically related) 
periodicities are detected for each RRAT. In addition, we find that pulse emission 
is a random (i.e. Poisson) process on short (hour-long) time scales but that most of 
the objects exhibit longer term (months-years) non-random behaviour. We find that 
PSRs J1819— 1458 and J1317— 5759 emit more doublets (two consecutive pulses) and 
triplets (three consecutive pulses) than is expected in random pulse distributions. No 
evidence for such an excess is found for the other RRATs. There are several different 
models for RRAT emission depending on both extrinsic and intrinsic factors which 
are consistent with these properties. 
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1 INTRODUCTION 

Rotating Radio Transients (RRATs) are neutron stars 
which were discovered o nly through their isolated pulses 
iJMcLaughlin et al.l 120061 ). Some, however, have later been 
detectable through periodicity searches. The average in- 
tervals between detected pulses range from a few min- 
utes to a few hours and pulses have durations between 

2 an d 30 ms. Thus far , ~50 RRATs have been identi- 
fied l|Hessels et all [20oi: iDeneva et all [2OO9I: iKeane et all 



2OIOI: iBurke-Spolaor fc Bailed I2OI0I : iKeane et al] I2OIII : 



Burke-Spolaor et all 201lh . including the original 11 from 



McLaughlin et al.l l|2006l ). Periods ranging from 0.1 to 
seconds have been measured for 29 of these sources. Pe- 
riod derivatives have been measured for 14, allowing in- 
ference of spin-down properties su ch as characteristic ages 
and surface dipole magnetic fields (JMcLaughlin et all l2009l : 
JLvne et al.ll2009l : JKeane et al.ll2011h . The periods and mag- 



netic fields of RRATs are larger than those of normal 
pulsars, but the distributions of other spin-down proper- 
ties such as spin-d own energy loss rate an d characteris- 
tic age are similar (jMcLaughlin et al.ll2009l ). Despite this 
overall trend, the properties of individual RRATs vary con- 
siderably. Four RRATs, including PSRs J1826-1419 and 
J1913-I-1330, have spin-down properties consistent with 
the bulk of the normal radio pulsar population and two 
others, PSRs J1317-5759 and J1444-6026, have proper- 
ties similar to normal, older pulsars. Four others, PSRs 
J1652-4406, J1707-4417, J1807-2557 and J1840-1419 
lie just above the radio 'death-line' fe.g. lChen fc RudermanI 
ll993l : IZhang et al.ll2007l V However, some have more unusual 
spin-down properties. PSRs J0847-4316, J1846-0257, and 
J1854-I-0306, lie in an empty region of P — P space between 
the normal radio pulsars and isolated neutron stars (XINS) 
and PSR J1819-1458 has a high magnetic field of 5 x 10^^ G. 
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Because of the difficulties in detecting these sporadic 
objects, the total Galactic population of RRATs likely out- 
numb ers that of normal radio pulsars iJMcLaughlin et alJ 
|2006| ). though it is possible that the pop ulations are evolu- 
tionarily related (JKeane fc Krameiil2008l ') . Several ideas have 
been presented about the nature of the emission from these 
objects. It could be similar to that responsib le for the 'gi- 
ant p ulses' observed from some pulsars (e.g. iKnight et al.l 
I2OO6I ). It could also be that the sporadic emission is related 
to t he fact that these objects are near the radio 'death-line' 
fe.g. lChen fc Rudermanll993l : [zhang et al.ll2007l') and/or ar e 
examples of extreme nulling fe.g. iRedman fc Rankinll2009 ). 
The phenomenon has also been a ttributed to the presence 
of a c ircumstellar asteroid belt (|Lill200q : ICordes fc ShannonI 
120081 ) o r a radiation belt as seen in planetary magneto- 
spheres (|Luo fc Melrosd 120071). Or, pe rhaps, some are tran- 
sient X-ray magnetars (e.g. IWoods e t al. 2005). Another 
idea is that their properties lie at the extreme end of the 
population of normal radio pulsars. Weltevrede et al. (2006) 
show that PSR B0656-I-14, a nearby middle-aged pulsar 
which emits pulses with energies many times its mean pulse 
energy, would be discovered as a RRAT source if it were 
farther away. RRATs may also be considered as an e xtreme 
case of mode changing (see, e.g., IWang et al.l l2007h where 
the on st ate is less than o r about one pulse period. Fur- 
thermore, |Lxne..et_aiJ (l201(J) have recently shown that many 
pulsars exhibit a two-state phenomenon in which varying 
pulse profile shapes are correlated with variations in spin- 
down rates and implied changes in magnetospheric particle 
density. These changes are quasi-periodic, with timescales 
ranging from one month to many years. It could be that the 
RRATs are similar two-state systems, in which the profile 
changes are so dramatic to make them undetectable in the 
more common state. 

Determining the time variability and/or periodicity of 
the RRAT pulses is therefore an important diagnostic of 
the RRAT emission mechanism. While the p ulse profile, 
and, i n most cases, pulse intensity changes of iLvne et al.l 
()2010 i) are quasi-periodic, the pulse intensity distributions 
of normal pulsars and giant -pulsing pulsars are believed to 
be random over time. On the other hand, nulling pulsars 
in general show on and off timescales of more than one 
consecutive puls e, indicating largely non- random distribu- 
tions (see, e.g., iRedman fc RankinI 12009 ) . Radio emitting 
neutron stars often show transient spin-down phenomena 
as well. For instance, glitches, or sudden increases in the 
spin frequency, have been observed from young pulsars and 
one RRAT ( PSR J1819-1458). One of the glitches from 
PSR J1819--1458 was accompanied by a 3.5ct increase in 
the pulse detection rate (jLvne et al.ll2009l ). Radiative events 
do not normally accompany the glitches of normal radio pu l- 
sars, but are quite common for magnetars (jDib et al.ll2008al ). 
This, along with the high magnetic field of PSR J1819-1458, 
hints at a relationship with magnetars and also provides ad- 
ditional motivation to examine the pulse rate variations with 
time for all RRATs. 

We search for periodicities and quantify the randomness 
of the detected RRAT pulses in several different ways. We 
first search for periodicities in the pulse arrival times on 
minutes-year long time scales and pulse detection rates on 
month-year long time scales using a Lomb-Scargle analysis. 
We then quantify the randomness of the RRAT pulse arrival 



times using Kolmogorov-Smirnov tests on seconds-year long 
time scales. The observations are described in Section 2, the 
methods and results in Section 3, and the conclusions and 
plans for future work in Section 4. 



2 OBSERVATIONS 

All eight sources discussed in this paper were discovered 
by McLaughlin et al. (2006) in a re-analysis of data from 
the Parkes Multi-beam Pulsar Survey (PMPS). We have 
ignored three of the original 1 1 RRATs as their pulse detec- 
tion rates are too low to perform this analysis. The discovery 
data were taken between Jan 1998 and Feb 2002 and follow- 
up observations began in Aug 2003 and are ongoing using 
the 64-m Parkes telescope. Most of the observations used 
the central beam of the multi-beam receiver with a central 
frequency of 1.4 GHz and a bandwidth of 256 MHz. A few 
observations used other frequencies; we ignore these in our 
analysis to ensure uniformity of pulse detection rates. The 
sources have been observed at between 27 and 89 epochs at 
1.4 GHz, with each observation 0.5—2 hr in duration (see 
Table 1). 

One important consideration in our analysis is the in- 
fiuence of the interstellar medium on our observed pulses. 
For all of these sources, the predicted diffrac tive scintillation 
bandwidths at 1.4 GHz are less than 1 MHz (jCordes fc Lazid 
[2003), making modulation due to diffractive scintillation 
unimportant. In Table 1, we list the predicted timescales 
for refractive scintillation at our observing frequency of 
1.4 GHz, estimated from the predic ted diffractiv e scintilla- 
tion timesc ales and bandwidths f rom Cordes fc L azio ( 200^) 
(see, e.g., iLorimer fc Kramen l|2005l )). These timescales 
range from 21 to 197 days. However, the actual timescales 
could vary significantly from those predicted. The pre- 
dicted mo dulation indices due to refractive scintillation 
(jLorimer fc Kramer 2005) range from 0.09 to 0.17, mean- 
ing these are expected to be relatively minor contributions 
to pulse rate variations. 



3 ANALYSIS AND RESULTS 

Pulse detection is performed by dedispersing the data at the 
dispersion measure (DM) of the RRAT and at a DM of zero. 
Then pulses are searched for in both time series above a 5(t 
threshold using the pulsar processing package SIGPROQj. 
Pulses which are brighter at the DM of the RRAT are likely 
to be from the source. We inspect the pulses visually by 
checking for pulse shape and pulse phase consistency to be 
certain of their astrophysical nature. For some epochs which 
have large amounts of radio frequency interference, we ap- 
plied the a bove procedure but with multiple trial DMs as de- 
scribed bv lMcLaughlin et al.l l|2009l) . If more than one pulse 
is detected within an observation, a second check based on 
the known period of the source can be made by requiring 
that all pulses have arrival times which differ by integer mul- 
tiples of the period. For the sources with phase-connected 
timing solutions, we check that the pulse arrival time is con- 
sistent with the solution. In Table 1, we list the number of 
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epochs for which pulses were detected for all sources. In Ta- 
ble 3 we list the total number of pulses detected within the 
entire time span of observations. 



3.1 Periodicity search 

The Lomb-Scargle test (Scargle et al. 1982) is a statistical 
procedure for uncovering periodic signals hidden in noise. 
We use this technique in our analysis as our data are un- 
evenly sampled, thereby making standard Fourier analysis 
di fficult. The implem entation we used (Numerical Recipes, 
see lPress et al.l (| 19861 )) utilises a version of the periodogram 
with modifications by Scargle (1982) and Home and Baliu- 
nas (1986). Applications of this test to radio pulsar data can 
be found in lBailes etlo] (|l997l ) and lKramer et all l|2006l ). 

The Lomb-Scargle test reveals signals in the power 
spectral density distribution of a source, with the presence 
of a sinusoid of certain frequency indicated by a peak in the 
spectrum at that particular frequency. The trial frequen- 
cies at which the periodogram is evaluated are chosen to be 
a finite evenly spaced set. For a time series X(ti) with A'o 
number of elements where i = 1,2,. . . TVo, the Scargle angular 
frequencies range from uj = 2n/T to u = ttNo/T (or peri- 
ods from T to 2T/No), where T is the total time interval. 
The searched frequencies therefore range up to the Nyquist 
frequency. The number of frequencies searched is obtained 
from the empirical formula (Home and Baliunas 1986) 



N, 



-6.362 + 1.193 X TVo + 0.00098 x iVo . 



(1) 



The likelihood of the existence of a signal or the level of 
significance is calculated as a detection threshold Zq (Scargle 
1982), 



Zo 



-ln[l-(l-po) 



i/Ni 



(2) 



The false alarm probability po is the probability that a peak 
of power Zo will occur in the absence of a periodic signal. 



3.1.1 Periodicity in pulse arrival times 

We have searched for periodicities in the pulse arrival times. 
In order to do this, we have created a time series by account- 
ing for all rotations of the pulsar during each observation 
and assigning a delta function (i.e one for each detection 
and zero for each non-detection). We performed this search 
on the entire time span of observations. 

In the time series for the entire time span of obser- 
vations, we searched over ~1, 000, 000 periods spaced as 
outlined at the beginning of section 3.1, ranging from 
~25 minutes to 2300 days. Six RRATs PSRs J0847-4316, 
J1317-5759, J1754-30, J1819-1458, J1826-1419, and 
J1913-I-1330, show significant periodicities in the arrival 
times on these timescales. The peaks of highest sig- 
nificance for PSRs J0847-4316, J1317-5759, J1754-30, 
J1826-1419, and J1913-f 1330 are at 3.8, 1.6, 1.4, 3.6 and 
11.3 hours respectively, while PSR J1819— 1458 shows a 
long term periodicity of 2102 days. There are many other 
significant periods and harmonics for these RRATs, with 
eight, three, one, thirteen, fifteen, and eight independent 
(i.e. non-harmonically related) periodicities with significance 
greater than 99% (2.5a) for PSRs J0847-4316, J1317-5759, 



J1754-30, J1819-1458, J1826-1419, and J1913-H330, re- 
spectively (See Table 2). These periodicities range from 
hours to years. Because of the large number of detected pe- 
riodicities, we do not list them all here. 

In order to determine the time dependence of the peri- 
odicities, we divided the time series to halves and quarters 
and performed the search again. For every RRAT, all pe- 
riodicities detected in the full series were re-detected in at 
least one quarter subsection with lower significances. None of 
the periodicities were re-detected in every quarter, though 
seven out of the eight non-harmonically related periodicities 
of PSR J1913-I-1330 were re-detected in three of the quar- 
ter datasets. Similarly five out of the 13 non-harmonically 
related periodicities of PSR J1819— 1458 were re-detected in 
three quarters with significances greater than 95%. The rest 
were re-detected in only one or two quadrants. 

All detectable (i.e. within the searched range) non- 
harmonically related periodicities of PSR J1913+1330 were 
re-detected in both halves of the dataset. For the rest of 
the RRATs only about five (for PSR J1819-1458) to two 
(for PSR J1913-I-1330) independent periodicities were re- 
detected in both halves. However every periodicity was re- 
detected in at least one half section of the dataset with lower 
significance. These results in general show that the periodic- 
ities persist throughout the entire time span of observations. 

In order to gauge the reality of the periodicities, we ran- 
domised the time series of detections and non-detections by 
placing the pulses randomly within the observation windows 
and repeating the analysis. We found no periodicities with 
significance greater than 30% in any of these randomised 
time series, which suggests the periodicities found are real. 
Figure 1 shows the power spectra for the pulse arrival times 
from the randomised time series for the eight RRATs. 

3.1.2 Periodicities in daily pulse detection rates 

We have also applied this method to look for periodicities 
in the daily pulse detection rates. For each day, the obser- 
vation length and the number of detected pulses were used 
to calculate the rate of pulse detection. Figure 2 shows how 
this rate varies for the eight RRATs. We then applied the 
Lomb-Scargle analysis to these rates, with the results of 
this analysis shown in Figure 3. We list the most significant 
period in Table 1 along with its significance. 

We have performed white noise simulations and Monte 
Carlo simulations to verify the significance of the period- 
icities. In white noise simulations, the daily rates were re- 
placed by random Gaussian noise. We could then calculate 
the power spectrum amplitude corresponding to the desired 
false alarm probability. In the second method Monte Carlo 
simulations were used to generate spectra from random time 
series which have the same sampling and the cumulative 
probability distribution of their maximum amplitude is cal- 
culated. We then fit this distribution to Equation 2, mini- 
mizing x^ to determine an effective value for Ni as this deter- 
mines the false alarm probability for a given power spectral 
density. These tests verified the significances that we have 
quoted. 

PSRs J1819-1458 and J1754-30 have periodicities 
with greater than 2a significance at 1260 and 994 days re- 
spectively, while PSRs J1846-0257 and J0847-4316 have 
periodicities with greater than la significance at 135 and 
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Figure 1. Power spectral density vs. period from the Lomb-Scargle analysis on pulse arrival times for PSRs J0847— 4316, J1317— 5759, 
J1444-6026, J1754-30, J1819-1458, J1826-1419, J1846-0257, and J1913+1330 (left panel) and the corresponding plots when the 
arrival times are randomised (right panel). The randomised time series do not show significant peaks in the spectra. The dashed line 
represents the 99% significance level. Note the different y-axis scales in the left and right panels. 
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Figure 2. Daily pulse detection rate vs. date for PSRs J0847-4316, J1317-5759, J1444-6026, J1754-30, J1819-1458, J1826-1419, 
J1846-0257, and J1913+1330. 



1040 days respectively. The remaining four RRATs do not 
show a periodicity of significance greater than the la level. 
The significances obtained for the peak power spectral den- 
sity for PSRs J1819-1458 and J1754-30 from white noise 
simulations are 97% and 96% respectively. The significances 
calculated from Monte Carlo simulations are 97% and 89%. 



These periodicities are different from the ones detected 
in their pulse arrival times except for the I26OI240 ^ay pe- 
riodicity of PSR J1819— 1458, which is close to the pulse 
arrival time periodicity of 1186±7 days. In general the pulse 
arrival time and pulse rate recurrence periodicities are ex- 
pected to be independent though it is possible for them to be 
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Figure 3. Power spectral density vs. period from the Lomb— Scargle analysis on daily pulse detection rates for PSRs J0847— 4316, 
J1317-5759, J1444-6026, J1754-30, J1819-1458, J1826-1419, J1846-0257, and J1913+1330. The dashed Une represents the 95% 
significance level. 



Table 1. Spin period, distance, time between the first and last epochs of observation, the number of epochs, fluctuation period of the pulse 
detection rate, the significance of the most pr ominent peak in the spe ctrum, predicted refractive scintillation timescale (using diffractive 
scintillation timescales and bandwidths from ICordes &: Lazid l|2002h '). predicted modulation index due to refractive scintillation, and 
base-10 logarithms of the derived parameters characteristic age and surface dipole magnetic field strength. 



Name Period Distance Data Span No Pi 

PSR (sec) (kpc) (days) (days) (days) 



AtRiss mRiss log[Tc] log[_B] 
(days) (yr) (G) 



J0847- 


-4316 


5.97 


3.4 


J1317- 


-5759 


2.64 


3.0 


J1444- 


-6026 


4.76 


5.5 


J1754- 


-30 


1.32 


2.2 


J1819- 


-1458 


4.26 


3.6 


J1826- 


-1419 


0.77 


3.2 


J1846- 


-0257 


4.47 


5.2 


J1913+1330 


0.92 


5.7 



2438 
2324 
2372 
2372 
2375 
2375 
2193 
1874 



69 


1040+^;^^ 


89% 


118 


0.11 


5.9 


13.4 


94 


560lf, 


44% 


31 


0.15 


6.5 


12.8 


72 


145 ±4 


5% 


197 


0.09 


6.6 


13.0 


48 


994lf,? 


95% 


21 


0.17 


- 


- 


72 


1260t^«g 


97% 


117 


0.11 


5.1 


13.6 


55 


162±5 


24% 


65 


0.13 


6.1 


12.4 


46 


135±4 


89% 


152 


0.10 


5.6 


13.4 


27 


178±9 


48% 


87 


0.13 


6.2 


12.4 



the same. In order to further gauge the rehabihty of our re- 
suhs for PSRs J1819-1458 and J1754-30, we created 1000 
random sequences by assigning the measured pulse detec- 
tion rates to randomly selected MJDs and calculated the 
number of times a peak of the same or higher significance 
appeared in the 1000 trials. This number was six and ten 
for PSRs J1819~-1458 and J1754-30, respectively, suggest- 
ing that the periodicities detected are real and that their 
significances may in fact be underestimated. However, the 
detected periodicities are at 40% and 50% of the total ob- 
servations lengths for PSRs J1754-30 and J1819-1458, re- 
spectively (see Table 1) ; longer observation spans are neces- 
sary to determine whether they are real. We also note that 
the significance of the peaks depends on the ranges of fre- 



quencies searched and that we have not searched frequencies 
higher than the Nyquist frequency or lower than 1/T. 

In order to determine whether there was any depen- 
dence of significance on period, we simulated sinusoidal sig- 
nals of various periods with additive random Gaussian noise 
and then applied the Lomb-Scargle algorithm. We found 
that the significance of the detected periodicities is inde- 
pendent of period. 



3.2 Randomness tests on pulse arrival times 

The Kolmogorov-Smirnov (KS) test (see, e.g.. Press et al. 
1986) is a statistical procedure which determines the degree 
to which two datasets differ. It compares the cumulative 
probability distributions of both datasets by calculating the 
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Table 2. Periodicity in the arrival time P2, and the number of harmonics (nh) detected for each periodicity in the entire time series 
for the six RRATs. Only the most significant harmonic for each period with a significance larger than 99% is given. The periodicities 
are listed in the order of their spectral powers. We required the ratio between integer multiples of periods to be less than 1.001 to be 
considered a harmonic. 



PSR J0847-4316 


PSR J1317-5759 


PSR J1754-30 


PSR J1819-1458 


PSR J1826-1419 


PSR J1913+1330 


P2 nh 


P2 nh 


P2 nh 


P2 nh 


P2 nh 


P2 nh 


(days) 


(days) 


(days) 


(days) 


(days) 


(days) 


0.163 10 


0.079 7 


0.066 2 


2101.98 1 


0.158 325 


0.479 751 


1.22 16 


1.16 1 




0.996 33 


0.188 311 


7.601 124 


265.74 1 


0.138 6 




1186.5 1 


0.869 143 


0.324 873 


3.08 16 


0.167 4 




1.22 34 


0.026 969 


1.624 437 


0.753 1 






0.907 38 


0.036 557 


0.618 674 


0.343 14 






0.776 39 


0.300 255 


34.657 28 


0.998 16 






0.548 36 


1.462 87 


0.244 938 


0.599 11 






4.668 26 
0.153 20 
3.050 13 
19.030 8 
6.942 9 
10.202 16 


1.239 100 
0.048 490 
0.407 213 
10.934 14 
7.353 20 
162.892 1 
327.001 1 
117.841 1 


1874.01 2 



maximum deviation 



max 

-oo<a;<oo 



\ C^{x) - C2{X) \, 



(3) 



of cumulative distribution functions (CDFs) C'l and C2. 
The accuracy of this test increases with the number of data 
poin ts and is expected to be accurate for four or more points 
([Stephens 197Q|). From this the probability that two arrays 
of data values are drawn from the same distribution can be 
calculated. Small values of this probability (i.e. large values 
of D) suggest that the distributions being tested differ. 

The test can therefore be used to explore whether the 
observed pulse sequences are consistent with random (i.e. 
Poisson) distributions. We have done this both by comparing 
the CDF of our data with the CDF of a simulated randomly 
distributed pulse sequence (i.e. kstwo test as implemented 
in Numerical Recipes) and by comparing the CDF of our 
data directly with the uniform CDF (i.e. ksone test). We 
have performed this test both for single days and for the 
entire time span of observations. For each single data set, we 
created simulated distributions by placing the same number 
of detected pulses at random times within the observation, 
with the constraint of allowing only one detected pulse per 
rotation. We also tested randomness on longer time spans 
by creating random time sequences of length equal to the 
total observation lengths. 

In Table 3, we list the total number of pulses observed 
from each RRAT along with the numbers of pulses detected 
on individual days. We also list the probabilities of the pulse 
sequences being random on a single day, averaged over all 
days of observation by comparing with simulated random 
pulse sequences and directly with a uniform CDF. The sim- 
ulated probabilities on each day are averages of 100,000 trials 
of different randomly generated data sets. We also list the 
rms deviations of the averages. Table 3 also lists the prob- 
ability of pulses being randomly distributed on long data 
spans for simulated data (again averaged over 100,000 re- 
alisations) and through comparison with a uniform CDF. 
Figure 4 illustrates this comparison for two RRATs. 

In general, the pulse distributions appear to be consis- 



tent with random (i.e Poisson) distributions. For simulated 
data, the single— day pulse distributions have probabilities 
ranging from 0.09 (for PSR J1754-30) to 0.52 (for PSR 
J0847— 4316) of being random when only days with more 
than two detections are considered. When days with more 
than three are considered the probabilities range from 0.16 
(for PSR J1826-1419) to 0.53 (for PSR J0847-4316), in- 
dicating that the bias for small numbers of pulses is small 
(there are no epochs with three or more pulses for PSR 
J1754— 30). When we compare directly with a uniform CDF, 
PSRs J1826-1419 and J1913-I-1330 show more evidence (i.e. 
probabilities of 0.05 and 0.06) for non-random pulse emis- 
sion. PSR J1754— 30 shows evidence for non-random behav- 
ior in both tests. As this is the RRAT with the lowest overall 
pulse rate, however, the tests are not very accurate. 

The total number of detected pulses from each source 
ranges from ~40 to ~1100, making our tests for randomness 
on long time scales very sensitive. When long time spans 
are considered, only PSRs J1444-6026 and J1846-0257 
are consistent with random distributions. The probabilities 
given by using the ksone and kstwo tests are similar. This 
is not surprising given the large numbers of pulses. We have 
tested these analysis techniques on intervals between pulses 
and checked for their consistency with an exponential distri- 
bution. Both pulse arrival time and pulse interval tests give 
identical results as expected. Furthermore we have searched 
for any correlation between the pulse intervals on different 
timescales and have found none. 

We explored the possibility of the telescope zenith angle 
at the time of the observation being responsible for varia- 
tions in pulse detection rates, as the system temperature 
depends on this factor and to account for effects based on 
varying gain and spillover. The KS test was carried out on 
pulse sequences for lower and higher angles (i.e. angles less 
and greater than 45 degrees) for this purpose. These results 
indicate that there is little dependence of the detected rates 
on the zenith angle, with probabilities ranging from 0.31 (for 
PSR J1317-5759) to 0.49 (for PSR J1819-1458). 

We have also searched for clustering of pulses. For all of 
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Table 3. Total observation time, total number of detected pulses, minimum, maximum, mean (with standard deviation in parentheses) 
of the number of pulses observed per day, average probability of short term randomness from the simulations (Pri), the rms of the 
average probability (cri), average probability of long term randomness from the simulations {Pr2), the rms of the average (a"2), average 
probability of short term randomness from direct comparison with the uniform CDF (Pui), and the average probability of long term 
randomness from direct comparison with the uniform CDF (Pu2)- The probabilities for the simulations were calculated for 100 trials. 
The duration of each observation was 30 minutes to 1 hour on average; as can be seen from N2 there was one very long observation for 
PSR J1819— 1458. The days with less than three pulses detected have not been included in P^i or the Pui probability calculation. The 
numbers outside and inside parentheses for P^i , cri and Pui indicate those for days with more than two/three pulses detected. 
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Figure 4. Cumulative distribution function for data over the entire time span of observations of PSR J1754— 30 (an example of a 
non-random pulse sequence) and PSR J1444— 6026 (an example of a random pulse sequence) compared with the CDF of a uniform 
distribution (dashed line). The figure also shows the point at which the KS statistic D is measured (vertical dashed— dotted line). Small 
values of D imply that the distributions are similar to each other. It can be seen that the first plot differs significantly from a uniform 
distribution, whereas the second docs not. The probabilities of randomness arc 0.001 and 0.31 for PSRs J1754— 30 and J1444— 6026 
respectively. 



the RRATs, the pulses usually occur singly with occasional 
consecutive pulses detected. We have measured the total 
number of doublets, triplets and quadruplets (i.e. instances 
of two, three, and four consecutive pulses) and compared 
with the number of doublets, triplets and quadruplets found 
in simulated random distributions. Table 4 lists the results of 
our analysis for 1000 simulated distributions. The number 
of doublets detected is higher than the number expected 



for a random distribution for both PSRs J1819— 1458 and 
J1317-5759. PSR J1819-1458 shows many more doublets, 
triplets and quadruplets than expected, with one instance 
of nine consecutive pulses. The mean duration of the on 
states (at which pulses are detected) for the eight RRATs 
in units of the period are also listed in Table 4. These were 
calculated by taking into account the number of detected 
multiplets and their duration. For the RRATs which do not 
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exhibit multiplets, we can only give an upper limit to the 
on-state duration of one period as the mean on-state will 
be less than one period. 



4 DISCUSSION 

We have shown that six of the RRATs have periodicities of 
significance greater than 99% in the pulse arrival times. The 
periods of the most significant peak range from 1.4 hours 
(for PSR J1754-30) to 2102 days (for PSR J1819-1458). 
No significant periodicities were detected upon randomiz- 
ing the time series, showing that these periodicities are real. 
We do not find any relationship between the number and 
significance of detected periodicities and spin-down proper- 
ties such as period or characteristic age. It is possible that 
some of the periodic behavior is due to refractive scintil- 
lation. However, the number and wide range of timescales 
of the periodicities are impossible to explain with refractive 
scintillation alone. 

The shorter timescale periodicities in pulse arrival times 
are similar to typically obser ved nulling timesc ales, which 
range from minutes to days (|Wang et al.l 120071 ). Explana- 
tions for pulsar nulling include an empty si ght-line pass- 
ing t hrough the sub-beam structure (Dcshpa nde fc Rankin 



200 ll ). a reversal of the emission direction (JMelikidze fc Gil 
20061 ). pulsar emission ceasing temporarily due to intermit- 



tent failure of pair pr oduction (Zhang et al. 20 07), an as- 
teroid belt of material dCord es & Shannon 200o) or changes 
in m agnetospheric currents (Lvno ct al. 2010; Wang ct al.| 
120071 ). Any combination of these could explain the extreme 
pulse-to-pulse variability of the RRATs and also the longer 
term periodicities. 

The significant periodicities found in the RRAT pulse 
arrival times may suggest a relationship with pulsars whose 
spin-down rates and pulse shapes un- dergo periodic varia- 
tions j_withjm^ed_changes in magnetospheric particle den- 
sity (jLvne et al.ll20ld ). The prototype of this source class 
is B1931-I-24 (Kramer e t al. 2006). The asteroid belt model 
of (jCordes fc Shannon! |2008 ) attributes the 40-day on/off 
timescale of this pulsar to an asteroid with eccentric 40-day 
orbit. It may be that the pulse-to-pulse variability of the 
RRATs is due to a similar process happening on very short 
timescales. However, because the on states of the RRATs 
are so short, it is impossible to measure period derivatives 
during the on and off states. It could also be that a simi- 
lar process is causing the longer term periodicities in pulse 
arrival times. 

If we apply the model of (jCordes fc Shannonl |2008) to 
the RRATs, multiple asteroids of an asteroid belt could be 
responsible for the observed periodicities in the arrival times. 
This is consistent with the large root mean square timing 
residuals which range from 1.1 ms (for PSR J1913-I-1330) 
to 11.2 ms (for PSR J0847-4316) of these RRATs as an 
earth -sized asteroid would ind uce residuals of the order of 
1 ms (jCordes fc Shannonll2008l ). However it is possible that 
much of these large residuals are due to pulse-to-pulse jitter, 
indicating that the true asteroid mass cannot be determined 
by the residuals only. 

The periodic fiuctuations in pulse arrival times could 
also be due to non-radial oscillations which drive differ- 
ent emission modes, as often seen in white dwarf stars 



I Rosen et al.l I2OIII ). The fundamental oscillation periods 
for neutron stars are expected to b e on the order of 
milliseconds fReiseneg ger fc GoldreichI Il992h to seconds 
( M cDcrmott et al. 19881 ) for g-modes. These are far too 
short to explain the multiple periodicities seen, but it is 
possible that we are observing the beat frequency between 
a non-radial oscillation period and t hat from a long er 
timescale process like those observed in lLvne et al.l l|2010l ). 

We also searched for periodicities in the daily pulse de- 
tection rates. Six of the RRATs do not show any signifi- 
cant periodicities in their daily pulse detection rates over 
timescales of months to years. The exceptions are PSR 
J1819— 1458, which exhibits a 1260 day period with a signif- 
icance of 97%, and PSR J1754-30, which exibits a 994 day 
period with significance of 95%. We detect a peak of higher 
significance only 0.6% and 1.0% of the time in 1000 trials 
in which the rates were randomly assigned to the MJDs 
for PSRs J1819-1458 and J1754-30 respectively, suggest- 
ing that the significance of the peak may be underestimated 
by the Lomb-Scargle algorithm. Given eight trials, if all of 
the RRATs had no periodicities, we would expect to find 
one periodicity with significance greater than 88%. How- 
ever, two RRATs with significances greater than 95% are 
not expected. The predicted timescales for refractive inter- 
stellar scintillation of 117 days for PSR J1819-1458 and 
21 days for PSR J1754— 30 are much smaller than the re- 
ported periodicities in pulse detection rates, indicating that 
these periodicities are not likely due to scintillation. How- 
ever, these periodicities are roughly half of the total data 
span. Therefore, a longer time span of observations is nec- 
essary to confirm them as significant. 

There is evidence for changes in period derivative as- 
sociated wit h changes in pu lse detection rate for PSR 
J1819-1458 (|Lvneet al.ll2009l ). The peak in the rate imme- 
diately follo wing t he first (and largest) glitch at MJD 53926 
(JLvne et al.l |2009|) hints at a correlation between glitches 
and emission properties. If confirmed in future glitches, 
this correlation might suggest a link with the magnetars, 
for which radiati ve changes often accompany glitches (e.g. 
iDib et al.l l2008bl ) or the class of mode changing pulsars 
whose period der ivative undergoes quasi-periodic changes 
(|Lvne et al.ll2010l ). 

All of the sources exhibit random pulse distributions on 
singl e days. T his is similar to that observed for giant pulses 
(Kru ght et al.l l2006). However, the pulses from the RRATs 
arc wider than those of observed giant p ulses and the pulse 
amplitude distributions are different (jMcLaughlin et al.l 
[2009; Mil ler et al.l l2010i). This therefore supports the idea 
of Weltrevrede et al. (2006), who suggest that the RRATs 
could be normal, but more distant, pulsars like B0656+14 
which emits very bright, narrow pulses and a distribution of 
weaker, broader pulses. 

Over longer time spans, the pulse sequences of six 
of the RRATs show evidence for non-random behaviour. 
These same RRATs show significant periodicities in arrival 
times. The external influ ence from, e.g., an asteroid belt 
ICordes fc Shannonl2008l ) could explain this through uneven 
distributions of material. However, it is difficult to relate this 
long term non-randomness to either normal or nulling pul- 
sars as, to our knowledge, this has not been explored for 
either of these source classes. We do not find any obvious 
relationship between the randomness and derived spin-down 
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Table 4. Measured number Nm of doublets, triplets, quadruplets vs. the expected number Ne, probability of occurrences under an 
assumed Poisson distribution Pp for eight RRATs and the mean duration of the on state in units of the period. The expected numbers 
are the means of 1000 trials of randomly generated fake distributions. 



Name 




Doublets 




Triplets 




Quadruplets 


On state duration 


PSR 


Ne 


N„, 


Pp 


N, 


N„> 


Pp 


No 


Nn, 


Pp 


(upper limit) 


J0847-4316 


1.06 


1 


0.35 


0.004 





0.13 








0.06 


1.007 


J1317-5759 


1.44 


7 


0.003 


0.009 


1 


0.161 








0.021 


1.036 


J1444-6026 


0.05 





0.949 








0.924 








0.901 


1 


J1754-30 


0.014 





0.981 








0.971 








0.961 


1 


J1819-1458 


35.42 


111 


2.5 X 10-** 


1.3 


15 


9.9 X 10-25 


0.04 


7 


3.8 X 10-45 


1.172 


J1826-1419 


0.06 


1 


0.028 


0.001 





0.957 








0.943 


1.016 


J1846-0257 


0.06 





0.902 








0.858 








0.815 


1 


J1913+1330 


0.97 


1 


0.306 


0.015 


1 


0.356 








0.358 


1.022 



properties such as age, period, and surface dipole magnetic 
field. 

Nulling pulsars tend to have on-off sta tes which persist 
over more than one pulse period (|Wang et al . 2007) whereas 
the large majority of the RRAT pulses occur singly, aside 
from a few pulses from PSRs J1819-1458 and J1317-5759. 
Redman & Rankin (2009) showed that the majority of pul- 
sars null non-randomly and several studies (e.g. Rankin & 
Wright 2008 and Herfindel & Rankin 2009) have revealed 
periodicities in the null cycles for some pulsars. The RRATs 
have quite different properties in these respects, but be- 
cause the nulls of some pulsars (e.g . B083 4-I-06, B1612-I-07, 
and B2315+21: lRedman fc RankinI (|2009l )) are random and 
single-pulse nulls are occasionally observed, we cannot ex- 
clude the possibility of the RRATs being an extreme case of 
nulling. This is supported by the similar ages and spin peri- 
ods of many RRATs and nulling pulsars and the strong cor - 
relation between nulling and pulsar age (|Wang et al.ll2007^ . 

In summary, we detect highly significant periodicities 
in the arrival times of six of the RRATs. We detect no 
highly significant periodicities in the long-term pulse de- 
tection rates for the RRATs, although we find tentative pe- 
riodicities for PSRs J1819-1458 and J1754-30. All of the 
RRATs show random behavior on a single day and most of 
the RRATs show non-random behavior on long timescales. 
Most of the RRATs emit pulses singly, but a few do show 
evidence for clustering of pulses. It is clear that there are 
periodicities in the pulse arrival times for these objects. The 
cause of these could be circumstellar material, non-radial os- 
cillations, or another process. Radio monitoring over longer 
time spans and observations at other wavelengths may be 
useful to further understand the reasons for the RRATs' 
unusual emission. More theoretical work and further studies 
of the randomness of pulse emission in normal pulsars with 
different properties is also necessary. 
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